©Two Mathematical Models of Nonlinear Vibrations 

Model parameters are fit to empirical vibration data. 
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Two innovative mathematical models 
of nonlinear vibrations, and methods of 
applying them, have been conceived as 
byproducts of an effort to develop a 
Kalman filter for highly precise estima- 
tion of bending motions of a large truss 
structure deployed in outer space from a 
space-shuttle payload bay. These models 
are also applicable to modeling and 
analysis of vibrations in other engineer- 
ing disciplines, on Earth as well as in 
outer space. 

The first model is denoted the ampli- 
tude-dependent stiffness (ADS) model to 
emphasize the difference between it and 
the classical linear harmonic-oscillator 
model, in which stiffness is a constant. The 
ADS model is embodied in the equation 

x + £x + K{x,x)x = 0, 

where * is the instantaneous amplitude 
of the oscillating position or modal co- 
ordinate, ^ is a damping parameter, 
and K{x,x ) is the ADS. 

In the initial outer-space application, 
the ADS was represented by the follow- 
ing nonlinear function: 

K{x,x) = a + bA(x,x) + cA(x,x) , 

where a, b, and c are constant parame- 
ters to be obtained by fitting the model 
to empirical amplitude-versus-fre- 
quency data, and A(x,x) is a modal am- 
plitude. The amplitude-versus-fre- 
quency data are obtained by means of a 
moving-window estimation technique 
in which one analyzes the instanta- 
neous vibration waveform during a 
time window of about 90 percent of the 


time-average vibration period. The am- 
plitude and frequency are taken to be 
those of a sinusoid that makes the least- 
squares best fit to the instantaneous am- 
plitude during the window (see figure). 
The window is then moved by about 2 
percent of the average period and an- 
other best-fit sinusoid is found. This 
process is repeated until a suitably repre- 
sentative sample of the vibration wave- 
form has been acquired. 

The modal amplitude is given by 


A(x,x) 


x 2 + 


K(x,x) 


where K(x,x) is any reasonable approx- 



A Nonlinear Decaying Waveform is approximated 
with a best-fit sinusoid during a moving window. 
The resulting sinusoidal amplitude and frequency 
data are collected from fits for the entire sequence 
of window positions and used to characterize the 
frequency versus amplitude of the nonlinear wave- 
form. These frequency versus amplitude data are 
then fit to an amplitude-dependent stiffness (ADS) 
representation. 


imation of K{x,x) . One can refine the 
approximation iteratively, starting from 
K(x,x) = a, then using the resulting 
value of A(x,x) in computing a value of 
A'(x,x) by use of the above equation for 
A'(x,x ) . 

The second model, denoted the mo- 
ment-expansion (ME) model, is embod- 
ied in the equation 

x + M(x,x) = 0, 

where the function M(x,x) is a moment 
expansion that captures damping and 
stiffness effects. The moment expansion 
is given by 
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M(x,x) = y y s P J p ij x t P, 

j= 0 *=o 

where both i and j range from 0 to 3, 
except that there is no (i,j) = (0,0) 
term. In the original outer-space appli- 
cation, the parameters py are obtained 
from (1) modal position and velocity 
estimates obtained from Kalman-filter 
states and (2) derived accelerations. 

In a test relevant to the original 
outer-space application, the ADS and 
ME models were compared with each 
other, with a linear model, and with a 
prior nonlinear model known as the 
Duffing model. The ADS model was 
found to yield the least error. 
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@ Simpler Adaptive Selection of Golomb Power-of-Two Codes 

The selected code-parameter value is within 1 of the 
optimum value. 
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An alternative method of adaptive se- 
lection of Golomb power-of-two (GP02) 
codes has been devised for use in effi- 
cient, lossless encoding of sequences of 
non-negative integers from discrete 
sources. The method is intended espe- 
cially for use in compression of digital 
image data. This method is somewhat 
suboptimal, but offers the advantage in 


that it involves significantly less compu- 
tation than does a prior method of adap- 
tive selection of optimum codes through 
“brute force” application of all code op- 
tions to every block of samples. 

A rather lengthy discussion of back- 
ground is necessary to give meaning to a 
brief summary of this innovation. For 
positive integer, m, the rath Golomb 


code defines a reversible, prefix-free 
mapping of non-negative integers to 
variable-length binary code words. 
Golomb codes are optimum for geomet- 
rically distributed sources (a model that 
frequently arises in image compression) : 
In the case of a geometrically distributed 
random variable, 8, the appropriately se- 
lected Golomb code minimizes the ex- 
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